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DYNAMICS OF A ROTATING CYLINDRICAL SPACE STATION 


SUMMARY 

Equations describing the reaction of a space vehicle 
to rotating machinery and to mass motion inside the 
vehicle are formulated. These equations include a 
component that has been neglected in most earlier 
studies , i.e., the angular impulse introduced by the 
moving mass. The dynamics of a cylindrical space 
station rotating around its major axis once every 10 
seconds are evaluated numerically. The long cylinder 
provides the large distance between the axis of ro- 
tation and the crew compartment needed to generate 
comfortable artificial gravity. But the rotating cylin- 
der is known to have little dynamic stability around 
its longitudinal axis. Nevertheless , the stability char- 
acteristics of this configuration are found to be 
acceptable provided there is as little as 0.2°/o asym- 
metry between the major and intermediate moments 
of inertia. Detailed investigations were performed on 
an analog computer. These show the effect of path 
and timing on the magnitude of the disturbance 
caused by a mass relocation. Also analyzed in some 
detail are ways to minimize roll , either by corrective 
motion of the astronauts or by a simple control sys- 
tem that uses the gyroscopic torque generated by 
rotating the body-fixed axis of a flywheel. 

SECTION I. INTRODUCTION 

The environment necessary for the life support of 
astronauts during extended space flights will prob- 
ably have to include artificial gravity. During periods 
of coasting, this will be generated by continuous 
rotation of the space station around an axis that is 
sufficiently remote from the crew compartment to 
exclude extreme variations in the centrifugal forces 
within the compartment. 

Our investigation will be limited to a station of 
cylindrical shape, since this is the geometry most 
likely to be used for early flights of long duration. 
As a further qualification, since the rotating cylinder 
is highly stable around its two axes of large moment 
of inertia, only the requirements for stability around 
the third axis, i.e., the longitudinal axis, will be 
examined. 


SECTION II. ANALYSIS 
A. The Undisturbed System 

In the absence of external torques, the angular im- 
pulse of the system remains constant. Let us assume 
that the angular impulse is aligned with the Z axis of 
the space-fixed reference coordinate system X, Y, Z 
(Figure 1). Using the transformation matrix from the 
reference coordinate system to the body-fixed co- 
ordinate system expressed in Euler angles <f>, 9 , we 
get the components of the angular impulse in the 
body-fixed coordinate system x, y, z (Figure 1):* 


2 . i 



«/r, 6, (f> as shown are positive angles. 

Figure 1 — Euler Angles, Transformation from Reference 
to Body-Fixed Coordinate System 


*THis approach eliminates the drift problems that occur when 
Euler's moment equation is used for the analog simulation 
of torque-free rotational motion. Cf. Whittaker, E. T., A Trea- 
tise on the Analytical Dynamics of Particles and Rigid Bodies , 
New York, Dover (1944), pg. 144. 
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or 

H x — — Hs\n9 

H y = H cos 9 sin </> (2) 

H z = H cos 9 cos <f> 

The angular impulse is related to the angular ve- 
locity by the vector equation 


We now express the angular velocities in terms of 
gimbal angles and rates: 

p = <j> — ip sin 9 . 

q=j=9 cos <f> + ip cos 0 sin (7) 

r = — 0 sin <f> + ip cos 9 cos <f> 


H= T- n (3) 

The inertia tensor / is given by a matrix 



If there is no mass motion inside the space station, 
the components of the tensor are constant in the 
body-fixed coordinate system. Calling the compo- 
nents of the angular velocity in the body-fixed co- 
ordinate system p, q, r and using Equations (2), the 
vector equation ( 3 ) can be written in components 

H x = hp — h y q — h z r = — H sin 9 

Hy = — hyp + i y q — iyzr = H cos 9 sin <j> (5) 

H z =~ hzP — lyzq + I zr=H cos 9 cos cj> 


Or, solving for <j>, 9, and ip, we get 
<j> = p sin 9 

9 = qcos<j> — rsin<£ ( 8 ) 

ip cos 9 = q sin <j> + r cos <j> 

Finally, introducing ( 6 ) into ( 8 ), we get 

D <£ = Di + (D 2 sin </> + D 3 cos <j>) tan 9 
D 9 = D 2 cos <p — D 3 sin <j> (9) 

D ip cos 9 — D 2 sin </> + D 3 cos <j> 

Equations (9) describe the attitude of a body with an 
angular impulse. This system of first-order differen- 
tial equations can be simplified if we choose a par- 
ticular body-fixed coordinate system that is aligned 
with the principal axes £, 77 , f The inertia tensor / can 
then be expressed by a diagonal matrix: 


Using determinants, Equations (5) are solved for p, 
q, and r 


where 



D = 


h - 

ixy 


hy hs\ 

iy hz , 


/ xz ly Z 


L 


Di = 


— H sin 9 hy hz\ 

H cos 9 sin <j> iy f vz 

Hcos9cos<f> — l yz l z 


Di = 


l x — Hsin 9 — hz\ 

— h y H cos 9 sin <f> — i yz 

— hz Hcos9coscj> h 


D s = 


h — ixy 

/ xy / y 

hz lyx 


H sin 9 

/-/cos#sin</> 

Hcos9cos<f> 


**Cf. Fifer, Stanley, Analogue Computation , Vol. IV, New York, 
McGraw-Hill (1961), pg. 1091. 



and equations (5) simplify to 

fi,, / £ = — H sin 0 

0^ = /-/cos0sin<£ (11) 

Qj c z=z H COS0COS<t> 

Rewriting ( 8 ) in our new notation, 

0 = f 2 7? cos <£ — H^sin^ 

4> = n i + isin0 (12) 

cos 0 = sin <£ -j- Cl c cos $ 

Finally, we get 

®= f(t-i) cosesin24> 



(13a) 
sin 0 
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or, using 

cos 2 $ = y -|- y cos 2 

1 1 

sin 2 = y — y cos 2 <t> 

we get 

*-Hi h+D+Hird™ 2 * 03b) 

i= T[-i + /!r + i + (r-c ) cos2 ®] sin0 

(13c) 

For rotational symmetry U 7) = I 1 ) 

0, i>, are constants 
and 

&£= constant 
O t? == A sin <t> 

Cl c — A cos 

These equations describe the nutational motion of a 
cylinder of perfect rotational symmetry. The length 
axis is moving on a circular cone with constant angu- 
lar and roll velocity — a well-known phenomenon. 

Let us now consider the case of a body with near 
rotational symmetry around the length axis and small 
angle 0: 

/ t -/„ = A/«l { 

0 <4- 

and 



Combining the first and the, derivative of the sec- 
ond of these equations, we eliminate 0 and get 

® = -t(- J— 9 ? 4 ' S| n 2 * 04 > 

/.e., the equation for the mathematical pendulum. 
For A/ positive the rotation around the £-axis is stable, 
and for Al negative the rotation around the 77-axis is 
stable, the period T increasing with decreasing A/. 


B. Effect of Mass Motion Within the System 

Next, how does a change in the mass distribution 
inside the station affect the motion of the station? 
The total angular impulse does not change because 
the mass motion is a result of internal forces. Nor 
does the center of mass move. But the origin of the 
body-fixed coordinate system does move and the 
inertia tensor changes. In addition, we must remem- 
ber that the angular momentum varies with location. 
A mass close to the center of rotation has less angular 
momentum than a mass located farther out. Thus a 
mass that moves away from the center of rotation is 
accelerated, thereby slowing the part of the station it 
moves to. (This is equivalent to saying that a moving 
mass in a rotating system experiences the Coriolis 
force, which causes a change in the angular velocity 
of the station if it acts at an arm. See Figure 2.) 



Figure 2 — Coriolis Forces Resulting from Mass Motion in 
Rotating Space Station 

Let us assume that a mass m moves from p o to p 
inside the station. Although the position vector of 
the center of mass in the space-fixed coordinate sys- 
tem does not change, the location R of the center of 
mass in the body-fixed coordinate system does 
move. Initially, we have 

Smipi + m po = 0 

With m moved to p the center of gravity has moved 
to R, given by 

M m R = 2 n ii pi -f- m p 

mi=£m 

= m p — mp 0 ( 15 ) 

where M is the total mass of the system. 
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The angular impulse H, which remains unchanged, 
is given by 

H = 2 m, Ti X Vi 

all m 

with Ti the position vector of m, and v t the velocity of 
nii (relative to the center of mass) measured in space- 
fixed coordinates. 

We write the equation for the angular impulse in 
body-fixed coordinates: 

5=jS — * 

yi— ^ - (n) =^Ti+TTxn 

=^jH — £ + nx(pi — R) 

H = ^ rru {pi — R) X (j& — R) - h 

2 m< (pi - R) X n x (pi - R) 

Using the formula a X (b X c) =5 (a*c) — c(a-b), 
we get 

H = 2m,(ft — R) X — 2mi(pi — R) X R 
+ 52mi(jR — R) 2 — 2m,- (jw — R)-[Q* (ft R)1 (16) 

The first two terms show the effect of the mass mo- 
tion; the last two terms, which depend only on the 
mass distribution, must correspond with Equation 
(3). Evaluating the first two terms further, we note^ 
that because ]Z = 0 for all m f - except the moving 
mass m, 

2m/(joi — R) X fi = m{p—R) X p 
and , according to the definition of R, 

2mi(pi — R)= 0 

Therefore the second term equals zero: 

2m f (j5i — R) XR = 0 

The equation for the angular impulse is now 

H = m(p — R) X/^H-^n (17) 

The x component of / * Q (see Equation 16) is given 
by 

(7 r iT) x =p 2m ( [(x,. — X)- + (y, — V) 2 +(z ( — Z) 2 ] 

— pSm; (x, —X) 2 

— q 2 m, (x, —X) (y f — V) 

— rjm, (x, — X) fz ( — Z) 

= p2m 1 [(y f — V) 2 + (z 4 — Z) 2 ] 

(x, —X) (y, — V) 

— r2m,(x < -X)(z 1 -Z) 


and the other components of /-Q have the corre- 
sponding form. We use these equations to evaluate 
the inertia tensor I: 

/x = 2m l [(y l -Y) 2 + (z ( -Z) 2 ] 

= 2m ( (y 4 2 + z ( 2 ) — 2V2m ( Y, 

— 2Z2m iZi + M (Y 2 + Z 2 ) 

= . - m (y„ 2 + z„ 2 ) + m(y 2 + z 2 ) 

~2Ym(y-y 0 ) 

-2Zm(z-z 0 ) + M(Y 2 +Z 2 ) 

and (see Equation 15) 

L=L, 0 + m (y 2 — y„ 2 ) (z 2 — z 0 2 J 

m2 m2 

-TT ( y-^ )2 -TT (2 - z » )2 (18a) 

where l x 0 is the x-component of / before the mass m 
is moved from p to p 0 . 

Since the body axes are the principal axes for p = p 0 , 
we get 

m2 

l xv = m (xy — x 0 yj — — (x — x 0 ) (y — y g ) (18b) 

M 

Similarly, we get 

ly =/„,» + m(x 2 — • x 0 2 ) + m(z 2 — z 2 ) 

~^-(x-x 0 ) 2 -^-(z-z 0 ) 2 (18c) 

M M 

I „ = h, o + m (x 2 — x 0 2 ) + m(y 2 —y 2 ) 

— (x — x » )2 — -77- (y — y 0 )2 dad) 

M M 

m 2 

Lz — m{xz — x„ z„) — — (x — x 0 ) (z — z 0 ) (1 8e) 

M 

ly* = m(yz — y 0 z 0 ) -!IL-(y — y J ( z — z ) (1 8f) 

M 

where I y 0 and \ z 0 are quantities corresponding to 

o' 

Equations (17) and (18) describe the dynamics of a 
rotating space station with no external torques. 

For our numerical evaluation of Equations (15), 
(17), and (18), we have used a model of a space sta- 
tion similar to one proposed in a study conducted by 
the Propulsion and Vehicle Engineering Laboratory, 
Marshall Space Flight Center, at the request of the 
MSFC Future Projects Office (Figure 3). The station 
is nearly cylindrical and rotates around its axis of 
maximum moment of inertia. The rate of rotation, Q, 
is 0.628 rad/sec and the crew compartment extends 
12 to 18 meters from the center of rotation. The re- 
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suiting artificial gravity, Q 2 *r, ranges from 0.48g to 
0.72g. Three astronauts weighing 90 kilograms each 
are free to move about in this area. The maximum 
moment of inertia, ! z , is 6,000,000 kg rrr. The mo- 
ment of inertia around the length axis, l u ., is 112,000 
kg m 2 . 

Several different mass 'motions were studied. Our 
study showed that a station with rotational symmetry 
around the length axis responds to mass disturbances 
with large roll motions. This is demonstrated on Fig- 
ure 4, which shows the results of a digital evaluation. 
A continuous roll motion is produced in response 
to a mass motion from x = 15m, 2 = 0 to x = 15m, 
z = 1.8m. The deflection in 0 is negligible. (The de- 
flection in 6 consists of a component proportional 
to cos and, according to Equation (13c), superim- 
posed on this a component proportional to </>. The 
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Figure 4a — <f> Response of Station 
to Mass Motion from x ~ 15m, z = 0 to x — 15m, z = 1.8m 
(Given Rotational Symmetry Around Length Axis) 


component proportional to cos <f> is caused by the 
angle between the principal axis $ and the body axis 
x.) The mass motion converts the z-axis into an axis 
of intermediate moment of inertia, the rotation 
around this axis is unstable, and because the motion 
(the distance of m from the center of rotation in- 
creases) introduces a small positive initial roll ve- 
locity, the roll motion is continuous. 

But even as little as 0.2% asymmetry is sufficient 
(for the model used in our numerical evaluation) to 
insure that the roll deflection in response to a single 
motion is limited to less than 10°. However, because 
of the absence of damping, whether or not a new 
mass motion will increase or decrease an existing 
oscillation depends on the timing of the new mass 
motion. (See the discussion of cancellation of roll 
by motion of the astronauts, page 34.) 



(Given Rotational Symmetry Around Length Axis) 
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Another disturbance that can be considered is the 
effect of rotating machinery (generators, pumps, etc.) 
on board the station. We modify Equation (17) to 
include the effect of rotating machinery with body- 
fixed axes: 

H=m(p — R) Xf + tt + h (19) 

where h is the angular impulse of the rotating 
machinery. Written in components, Equation (19) 
becomes 

— Hs\n6 = h x + pl x —ql xv — rl xz 

+ mUy-Y)z-(z-Z)y] 

H cos 6 sin <£ = h y — p I XI) + q /„ — r t yz 

+ m [(z — Z) x — (x — X)z] 

H cos 6 cos$~h. —p l xz — q l yz + rl z 

+ m [(x — X) y — (y — Y)x] 

We see from these equations that flywheels can en- 
force the correct attitude in the steady state, in spite 
of mass displacements. 

In the steady state, with p, q, x, y, and z = 0, <f> and 
6 can be reduced to zero if 

h x = rl* z 

h „= rl V z 

We can determine the required angular impulses to 
do this for the maximum possible values of l xz and 

v 

l xz =mxz = 270kg • 18m • 1.8m = 8,750kg m 2 

l yg = myz = 270kg • 1.27m ♦ 1.27m = 435kg m 2 

Because the roll deflection is considerably larger 
than the pitch deflection, it is more important to re- 
duce roll. Fortunately, this requires a twenty times 
smaller angular impulse. Moreover, because a de- 
flection in the roll orientation is less disturbing to 
the astronauts than a roll oscillation, only methods of 
limiting roll oscillation were studied. Two methods 
were studied in detail: controlled motion of the 
astronauts and automatic control by flywheels. 

To simplify our analysis of the effect of rotating 
machinery with body-fixed axes, we shall assume 
that the body-fixed axes are the principal axes and 
that the astronauts are at rest. Equations (19) then 
become 


— H sin 6 = h x + pl x 
H cos 6 sin </> = h 9 -\-ql v 
H cos 0 cos cj> = h z + r l z 


We can estimate the initial roll velocity caused by 
a fast start (assuming that 6 = 0 initially) from 

t x P + h x = 0= h x -\- l x $ 0 

If we limit the roll deflection to a maximum of 0.1 
rad, the roll oscillation will be approximately sinu- 
soidal, with the relation between the maximum 
angular velocity and the angular deflection given by 


It follows that 

| h x | < izZlL -o.l rad = 1 850kg m 2 sec- 1 

This angular impulse causes a deflection in 6 after 
the roll oscillation is damped out: 

q = _K 1850 

H 6,000,000 • 0.628 

= — 0.5 • 10“ 3 rad 

But, because of their small size, 6 and 0 can be neg- 
lected, and we get for the components of O: 

q = ip cos 6 sin <j> r = ^ cos 6 cos 

Using these expressions, we get 

h z = (H — ij/ l z ) cos 6 cos <j> 
or 


Thus, because H is very large, the change in ^ re- 
sulting from rotating machinery with its axis parallel 
to z can be neglected. 

Finally, for h y , we get 

h y ~ H ^ 1 — cos 6 sin <j> 

Hence an angular impulse h y causes a steady-state 
deflection in <j>. Imposing the same limit on the 
steady-state deflection as we did on the roll oscil- 
lation, we get 


h v< H (l — Y)' 01 
< 6,000,000 • 0.002 • 0.628-0.1 
< 750 kg m 2 sec- 1 


The magnitude of h x that can be tolerated is de- 
termined by the transients it introduces into roll 
during starting and stopping. 


In sum, then, rotating machinery interferes least 
with the dynamics of the rotating station when its 
axis is aligned with the z-axis of the station. 
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SECTION III. NUMERICAL EVALUATION 

A. The Analog Computer Program 

To study the behavior of Equation (19), we pro- 
grammed it on an analog computer. Mechanization 
in the most general manner had to include the fol- 
lowing: the variation in the nine components of the 
inertia tensor resulting from changes in the mass 
distribution, the angular impulse of the moving mass 
inside the station, and the relation between the com- 
ponents of the angular impulse (in the body-fixed 
coordinate system) and the Euler angle rates. Seven 
multipliers, with 35 products, and two resolvers were 
required. 

By various linearizations and approximations, this 
requirement was reduced to four multipliers, with 11 
products. The largest error (7% in the worst case) 
was introduced into the factor of </> by the replace- 
ment of 1 — ljl z with 1 — l ytO /l gf0 - 


The following simplifying assumptions were made: 
that m << M, that the position of m in the undis- 
turbed system is at y Qt z„ = 0, that changes in l x , l yf l s 
resulting from mass displacements can be neglected, 
and that small angle approximations will have no 
significant effect on the results. We therefore get 

l xy ~ rnxy; f yz zz myz; i zx 7Z mxz 

• • * • * f-f 

— ^ 0; q ~ 0 + ^ <£; r ~ ^ ~ — 

' z 

l x <}>~ — H‘(l — y-) 6 + mxy0 + mj-xy^ 

H 

+ m--xz+m(yz — zy) — h x 

' z 

l v 6 zz H f 1 cf> + mxy — m —xy 0 -{- m — yz 

+ m (xz — zx) — h y (20) 

The analog computer program is shown in Figure 5. 


o o o o o o 
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Figure 5 — Basic Analog Computer Program 
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SPACE STATION THREE ASTRONAUTS 


The effect of a mass displacement depends on the 
path and on the schedule of the displacement. Fig- 
ure 6 shows the response of the station to the mo- 
tion of three astronauts (3 -90 kg = 270kg) from 
(12m, 0, 0) to (18m, 0, 1.8m) along three different 
paths: 

Path I. From (1 2m, 0, 0) to (18m, 0, 0) and then to 
(18m, 0, 1.8m) 

Path II. From (12m, 0, 0) to (18m, 0, 1.8m) and then 
to (18m, 0, 1.8m) 

Path III. From (12m, 0, 0) to (12m, 0, 1.8m) and then 
to (18m, 0, 1.8m) 

Check runs on a digital computer show good agree- 
ment, indicating that the approximations used have 
no significant effect on the results. It follows from 
Equation (13a) that the pitch amplitude is 0.007 of 
the roll amplitude. This too is in agreement with the 
computer results. The steady-state pitch deflection 
can be determined from Equation (E) of the Appen- 
dix, where the orientation of the principal axes of 
the space station relative to the body-fixed axes is 
described. 

PATH I 


B. Effects of Astronauts' Motion 

The walking speed was assumed to be 0.9 m/sec. 
To simulate changes in velocity realistically, a first- 
order lag with a time constant of one second was 
added to the velocity term. 

Motion in the yz plane, for example the motion 
of three astronauts from (18m, 0, 0) to (18m, 1.2m, 
1.2m), changes the direction of the principal axes 77, 
£ relative to the body axes and causes a roll oscilla- 
tion, but the deflection in this case is only 2°. If this 
transfer occurs in two phases, for example from (18m, 
0, 0) to (18m, 0, 1.2m) and then to (18m, 1.2m, 1.2m), 
a small initial roll velocity (+2.5 * 10 -3 rad/sec) is 
introduced in the second phase. The same two-phase 
transfer in a different sequence, /. e., from (18m, 0, 
0) to (18m, 1.2m, 0) and then to (18m, 1.2m, 1.2m), 
produces a small initial roll velocity of the opposite 
sign. 

Because starting and stopping occur in time inter- 
vals that are small compared with the period of the 
oscillation, the effect of these on the roll amplitude 
can be neglected (Figure 7a-c). 

PATH II PATH III 



Figure 6 — Response of Station to Motion in Unison of Three Astronauts from (12m, 0, 0) to 18m, 0, 1.8m) along Three 
Different Paths. (Note 100-to-One Scale Change in $ between Path I and Paths I! and III.) 


8 






Figure 7 — Response of Station to Motion in Unison of Three Astronauts (a-c) in yz Plane 
(x = 18m), (c, d) in xy Plane (z = 0), (e-g) Parallel to but Not Along x Axis, (h, i) 
in Circle in yz Plane (x = 15m) 


Motion in the xy plane (z = 0) causes no roll, as 
can be seen in Figure 7c, d. The ladder for going from 
a lower to an upper deck of the space station should 
therefore be located in the xy plane (z — 0). 

Motion parallel to, but not along, the x-axis (con- 
stant Zy^O) introduces large roll disturbances, mainly 
because of the change in angular momentum with 
distance from the center of rotation. The transfer of 
this momentum to the roll axis depends only on 
z, /. e., the torque generated by the Coriolis force 
(Figure 7e-g). 

Walking in a circle in the yz plane (Figure 7h, i) 
introduces two disturbances: the attitude of the prin- 
cipal axes relative to the body axes is varied, and an 
angular roll impulse is introduced. The first disturb? 
ance has a period corresponding to the time it takes 
to walk around half the circumference of the cylin- 


drical station. Its effect can be seen in <£. The second 
disturbance is the reaction of the station to the ini- 
tiation of circular motion by the astronauts. The 
initial roll velocity that results is given by 

l x $ o = mp*v— 270kg * 1.8m * 0.9 m/sec 

The amplitude of the resultant roll (7= the period 
of the roll motion) is 

«£= = 24-10- 3 rad 

v 2 TT 

All these disturbances are moderate. The roll 
angle never exceeds 10°. However, whether the roll 
of the space station is increased or decreased 
depends on the phase of the recurrence of a 
disturbance. 
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Figure 8 — Effect of One Astronaut Moving in x Direction Synchronously with Roll: (a) Check 
Run, Repeat of Figure 6, (b) Roll Unaffected by Motion in xy Plane (z = 0), (c) Damp- 
ing of roll (z > 0), (d) Augmentation of Roll (2 < 0), (e, f) Rofl Unaffected by Motion 
in Phase with 


Figure 8 shows the effect of one astronaut moving 
in the x direction synchronously with a roll oscilla- 
tion. Whether the effect is one of buildup or of de- 
cay depends on the phase relation. Hence it should 
be feasible to use the astronauts for corrective action 
to counteract roll of the space station. Figure 8b 
shows again (as did Figure 7c, d) that motion in the 
xy plane (z = 0) does not affect roll. Figure 8c shows 
a damping effect (z > 0); Figure 8d shows an aug- 
mentation (z < 0) for motion in the x direction in 
phase with </>; Figure 8e, / shows that motion in the 
x direction in phase with </> has no effect. 

A more effective and convenient way to reduce a 
roll oscillation is shown on Figure 9. Here the motion 
is parallel to, but not along, the y-axis (constantz 7 ^ 0 ) 
or parallel to, but not along, the z-axis (constant 
y^0) in phase with the roll velocity.* The same mo- 
tion in the opposite phase causes a buildup. 


Very little has to be added to the basic analog 
program to simulate motion of an astronaut syn- 
chronized with a roll oscillation of the station. The 
additions to the basic program are shown on Figure 
10. With Integrator 8 (Figure 10) providing an adjust- 
able first-order time lag, Amplifier 9 generates the 
sign of <£ or, with 51 closed, the sign of p. The bridge 
limiter inserts a threshold, so that corrective action 
is taken by the astronaut only if roll exceeds a certain 
limit. The x c , y r , and z r , commands are connected to 
their integrators over a pair of relay switches that 
keep the resultant motion within preset limits.* For 
example, the astronaut cannot move beyond the wall 
of the soace station, /. e., y and z must always be less 
than 1.8m. 

*Measured as <£, /. e., as the gimbal rate of the inertial platform, or 
as p by a rate gyro mounted on the frame of the space station 
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Figure 9 — Correction of Initial Roll Disturbance by Controlled Motion of One Astronaut 
in yz Plane 




Figure 10 — Computer Program: Motion of Astronaut Synchronized with Roll Oscillation of Station 




FLYWHEEL SPACE STATION THREE ASTRONAUTS 


C. Automatic Control of Station's Roll Motion 

As already mentioned, properly controlled varia- 
tions in flywheel rates could be used to damp out 
roll oscillations. We shall first consider the effect of 
a flywheel with its axis parallel to the x-axis. 

For simplicity, we assume that f xy 
that the astronauts are not moving. The station has 
a roll oscillation introduced by some previous dis- 
turbance. We use small angle approximations and 
Equation (19) becomes, with l z — l y = A* and /* > > 

Ah 


lyz, /«= 0 and 


Hd^h x + pl x 
H<j> ~ ql v 
H^rL 


r ~ 'A 


Thus 


and 


— — 

V /, ij L 


, A',, 

' 17 * 




H- r 2 J, _ hjL 
L L 


( 21 ) 


This equation shows that only a rate of the angular 
impulse of the flywheel, h xf generates a torque. Thus 
if fi x can be varied synchronously with <£, the roll 
oscillation of the station can be damped. But h x can- 
not be used to counteract a steady-state deflection 
in <j> caused by ! yz . However, as Figure 7a shows, the 
largest possible steady-state deflection due to l ye is 
only 2°. 



+ 18 — 

x (m) 

0 — 

y (m) 

+15 — 

-1 5 — 

z (m) 

+ 15 — 
0 — 
1.5 — 


4-25 I0 3 — 

6 (rad) 

o 

</> (rad) 

- 25 I0 3 

+ 25 — 

0 — ■ 


- 25 — 
+ 025 — 

17) 

0 — 


- 025 

— 

+ 200 — 


h c (kg m 2 sec -1 ) 0 


-200 - 
+ 200 - 


h y (kg m 2 sec -1 ) 0 — j 


!r 


i r 


L 


I 2 = 6,000,000 kg m 2 h c = 130 kg m 2 sec" 1 
I x = 1 12,000 kg m 2 


T_r 





VW*i /WV- 




— \AAAj — VW 1 — /W — — v - — v v* 


n i i i i i i i i i i i r i i f r 

100 200 300 400 500 600 700 800 900 1000 MOO 1200 1300 1400 1500 1600 1700 

TIME (SEC) -» 

Figure 11 — Correction of Initial Roll Disturbance by Flywheel Control System 
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Figure 12 — Computer Program: Synchronous Motion of Flywheel for Automatic Roll Control 


Another possibility is to use a flywheel with’ its 
axis parallel to the y-axis. Using the same approxi- 
mations, we get 

^^ — AL r 2 <£ + r^ (22) 

^ X f X 

We now get a damped roll oscillation if the angular 
impulse of the flywheel, h y/ is varied synchronously 
with — <£. Comparing Equations (21) and (22), we 
can see that for an equal damping effect 



An effective schedule for h was found to be a 
linear increase with time followed by an instanta- 
neous braking action (Figure 11). For this schedule 
hr is related to h x as follows (7 = the roll period): 



0.628h 


^max 


= 12h. 


'max 


The simulation assumes that the flywheel* is driven 
by a motor with a combined time constant of 10 sec- 
onds. (Therefore h y does not remain constant, as was 
assumed above, but decays during the cycle.) The 
voltage applied to the motor is constant; only the 
polarity is controlled. Brakes are applied before the 
voltage is reversed. Thus 

f = const, sign <f> = c * h c 

h = Jh — 

* 1+IOs 

The maximum required power, P, is half the maxi- 
mum torque multiplied by half the maximum angular 
rate, e.. 


With l F — 2kg m 2 and h c = 100 lb ft sec” 1 
= 130kg m 2 sec -1 , we get 


Control by varying h y thus turns out to be consid- 
erably more effective. 

The additions to the basic analog program to simu- 
late such a control system are shown on Figure 12. 
The discharge of Integrator 3 (Figure 12) through the 
two relay switches occurs in the intervals when the 
angular impulse of the flywheel, h c on Figure 12, 
approaches zero (adjusted by potentiometers 51 
and 52). 

The performance of the flywheel control system 
is illustrated in Figure 11. No phase adjustments were 
made to optimize the control. 


1 130 130 

max ~ 4 10 * 2 

= 211 watts 

The maximum angular rate is therefore 

w= =65 rad/sec 
2 

= 620 rpm 

*A hydraulic system with circulating pump and valves could well 
be another feasible solution to the control system problem. 
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Or, rather than building up an angular impulse 
and then wasting the stored energy by braking, we 
can use a gimbal-mounted gyro with its axis aligned 
initially with the 2 -axis (Figure 13). The required 
angular impulse, h y , is generated by rotating the gim- 
bal (its axis is parallel to the x-axis) through an 
angle 8: 

h |f = — / # O a sin 8 

(l 8 is the angular impulse of the spinning gyro.) 



The other two components of the angular impulse 
generated by this rotation are 

k=>g,J 

h z='» COsS - 

(Iq, x is the combined moment of inertia of the gim- 
bal and the flywheel around the x-axis.) Because it 
it small, the influence of these other two components 
on the motion of the station can be neglected. 

The torque generated by turning the gimbal is 

r=/7+n X h 

= Vo, x + *A (q cos 8 + r sin 8)] * / 

— [I 8 n 8 h cos 8 + p/A cos 8 + rl GtX h] ■ / 

— [/A* s ' n ^ “f - P^A s ' n s + * k 

The / and k components of L, along with the centrifu- 
gal forces, constitute the load on the gimbal — and 
spin — bearings. Sufficient power must be provided 
to replace the energy dissipated by friction. 

Because the station is rotating, L x =£0. Power is 
therefore required to turn the gimbal: 

p=i m 'h 


This power builds up potential energy: 

6 — 0 6 

• r-i .i r 

E = fL x S dt X , 8 2 J + /. • O, J r sin 8 d S 

6 = 0 0 

(q is neglected as small) or, considering r to be in- 
dependent of 8: 

£ l a ' 0/(1 — cos 8) 

and for 8 = —, 

2 

E ^ 82 watt sec 

This energy has only to be built up. As soon as the 
torque is removed, the gimbal swings back to the 
opposite deflection. 

It may be advantageous to use moderate rates of 
rotation of the gyro. The stored kinetic energy would 
thus be small, allowing the control to be switched 
on and the gyro to be started as required. The stored 
kinetic energy for an angular impulse, l 8 * Q 8 , of 
130kg m 2 sec -1 is, given that / g = 2kg m 2 , 

£ == 4250 watt sec 

SECTION IV. EFFECT OF EXTERNAL TORQUES 

Still to be mentioned are the effects of external 
torques. These have not been considered in our 
analyses mainly because external torques in space 
are small (gravity gradient, magnetic fields, radiation 
pressure, drag). They therefore change the angular 
impulse only very slowly: 



The same can be said of the impact of meteoroids. 
Roll motion resulting from such impacts can be 
damped out as it occurs and the axis of rotation will 
be in line with the direction of the angular impulse. 

Finally, it may occasionally be necessary to use jets 
to restore and maintain the desired rate of rotation 
and the direction in space of the axis of rotation. 

SECTION V. CONCLUSIONS 

Mass motion in a rotating space station of cylindrical 
shape causes the station to oscillate in the roll direc- 
tion. The more nearly symmetrical the station is 
around its length axis, the larger the amplitude of 
the roll oscillation will be. But, as is shown in the 
example investigated, a difference even as small as 
0.2% between the moments of inertia of the major 
and the intermediate axes suffices to keep the oscil- 
lation small. In addition, excessive roll can be re- 
duced by control systems using a flywheel or correc- 
tive motion of the astronauts. 
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APPENDIX — Orientation of the Principal Axes of the 
Space Station Relative to the Body-Fixed Axes 
The angular Impulse Is related to the angular ve- 
locity by 



-L, 


-I 


-/ 


-I, 



(A) 


For H pointing in the direction of the principal axis. 


tf = AO 


Writing thisJn components, 

/,P — /*,< 1 — I„r = *-P 

— i xi/ p + i„q — i u * r = kc i (B) 

-'«P-^q + ^=Ar 
It follows that the determinant 


~Ly Vy-H -l„ 


=0 


This cubic equation in A has three roots, A,, A 2 , A 3 , 
which correspond to 1^, l 1f/ l r The value for A closest 
to l x is l $ . Using A 1 = / £ , Equations (B) become 


— p — Uq — /« r =° 

Ly P + (fy /f) 9 lyz r=0 

LzP lyzP + ( h — 1 *) r= 0 
The first equation can be interpreted as a scalar 
product a • £2 = 0 indicating that a* is perpendicular 
to £2 with 

3 = (',■ — >*) ' ~ LJ — f x Jk 

The second equation yields a vector, 6, with 

b = ~ l ,y i + Vy — IJi — lyzk 

The vector aXb has the direction of the principal 
axis f. 

! i i k 1 
axb =!(/,-/,) , 

I 

We can now express the orientation of the princi- 
pal axes relative to the body-fixed axes by determin- 
ing the three Euler angles «, ft, y. We use the same 
sequence of rotation as we used in converting from 
the reference axes to the body-fixed axes. The z-axis 
agrees with y, the £-axis with «. We can therefore use 
the matrix of Equation (1): 


T cos /? cos y 

£2 7; = cos y sin /3 sin « — sin y cos a 
.£2 f J Lcos a sin /? cos y + sin a sin y 

P 
<7 
r 


cos ^ sin y — sin £ 

sin a sin /? sin y + cos a cos y sin a cos ft 
cos a sin ^ sin y — sin a cos y cos a cos /? 


(C) 


Because a X b has the direction of the £-axis, its com- 
ponent /*aXb is proportional to cos y cos ft and its 
component /-aXb is proportional to sin y cos /?. We 
determine y for the £-axis (selecting the smallest y) 
from 


/aXfa = sin y coSjS =tan /,»(/, — /*) + l„ • > xv 
iaXb cos y cos (S t yz l xy + l xz (t y — l t ) 

(D) 

We determine /? from 

**><£ = - si - n -.£_ = _ J-tan/3 (E) 

/aXb COS y cos ^ COSy 


u /„ + /„(/,-/«> 

To determine a, we rewrite Equations (2) using 
A 2 = It), 

(/. — h)p — l xv q — l xz r = o = c-n_ 
— UP - f (/, — l v )q — l vz r = 0 = d n 


and we use the relationship 




COS 

COS 


— tan « 
0 


Because /?, y are small. 


/cXd 

kcXd 


tan « 


(F) 
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